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On the Well-Posedness of the Vacuum Einstein's 

Equations 

Lavi Karp * 

Abstract 



p ! . The Cauchy problem of the vacuum Einstein's equations aims to find a semi- 

metric g Q( g of a spacetime with vanishing Ricci curvature R Qj/ g and prescribed initial 
y-n ' data. Under the harmonic gauge condition, the equations R a> ^ = are transferred 

into a system of quasi-linear wave equations which are called the reduced Einstein 
equations. The initial data for Einstein's equations are a proper Riemannian metric 
h a b and a second fundamental form K a 6. A necessary condition for the reduced 
Einstein equation to satisfy the vacuum equations is that the initial data satisfy 
Einstein constraint equations. Hence the data (h a b,K a b) cannot serve as initial data 
00 \ for the reduced Einstein equations. 

Previous results in the case of asymptotically flat spacetimes provide a solution 
■ to the constraint equations in one type of Sobolev spaces, while initial data for the 

| evolution equations belong to a different type of Sobolev spaces. The goal of our 

' work is to resolve this incompatibility and to show that under the harmonic gauge 

^ . the vacuum Einstein equations are well-posed in one type of Sobolev spaces. 

> 

X ' 1 Introduction 

This paper deals with well-posedness of the Cauchy problem of Einstein vacuum equations 

R Q/3 (g) = 0, a, (3 = 0,1, 2, 3. (1.1) 

Here R aj g(g) denotes the Ricci curvature tensors of a Lorentzian metric g. The unknowns 
are the coefficients g a/ g of the semi-metric g. 

The initial data consist of the triple (M, h a b,K a b), where M is a space-like manifold, h a b 
is a proper Riemannian metric on M and K a b is its second fundamental form (extrinsic 
curvature). The semi- metric g a p takes the following data on M: 

gOO | AT — — 1) gOa| A/ = 0, goi Af = h a & 



1^7 _ ^ bO01M ao a, b =1,2,3. (1.2) 
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The rest of the data dog Q o are determined through the conditions 




(1.3) 



Here denotes the Christoffel symbols and g^ 7 the inverse matrix of g^ 7 . Condition 
(1.3) is known as the harmonic gauge. Since (1.1) is a characteristic (see e.g. [14]), it is 
impossible to solve it in the present form. However, under the harmonic gauge (1.3), the 
vacuum Einstein equations (1.3) is equivalent to the reduced Einstein equations: 



Since g a( g has a Lorentzian signature, this is a system of semi-linear wave equations. The 
expressions Q 7 5(g,(9g) are quadratic functions of the semi-metric g Q/ g and its first order 
partial derivatives dgap. 

It is well known that the initial data (h a b,K a fc) cannot be prescribed arbitrarily and, in 
fact, the Codazzi equations, which relate the curvature of the manifold (V, g a /?) to the one 
of (M, h a b), lead to the Einstein constraint equations: 



Here R(h) and D a are the scalar and the covariant derivative with respect to the metric 
Kb- 

The fulfillment of the constraint equations (1.5) is a necessary condition for the solution of 
the wave equations (1-4) with the initial data (1.2) to satisfy the vacuum Einstein vacuum 
equations (1.1). We refer to [1], [24], [16] or [27] for a discussion of this fact. 

Conclusion: We conclude that any solution of the Cauchy problem for the Einstein 
equations (1.1) includes the treatment of the following two problems: 

(i) Solutions to constraints (1.5), which can be reduced to elliptic equations; 

(ii) Solutions to the reduced Einstein equations (1.4) with the initial data (1.2). 

We will deal with these two problems in asymptotically flat manifolds (AF). A Riemanian 
3-manifold (M, h a &) is (AF) if there is a compact subset K such that M\K is diffeomorphic 
to IR 3 \ -Bi(O) and the metric h a b tends to the identity e a & at infinity. 

1.1 Existence of the evolution equations 

The unknowns g aj g are functions of t = x° and x a , a = 1, 2, 3. We assume that for t = 0, 
the initial data (1.2) are given on a space-like (AF) hypersurface M ~ M 3 and denote the 
Bessel potential space on M 3 by H s . 



g af} d a d(3g 7 s = Q 7 5(g,<9g). 



(1.4) 




(Hamilton constraint) 
(momentum constraint) 



(1.5) 
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The following classical result was first proved by Y. Choquet-Bruhat [7] for s > 3 and 
improved by T. Hughes, T. Kato and J. Marsden [17]. 

Theorem A. Assume the following hold: 

(g Q/3 (0) - m Q/3 ) e H S+1 (R 3 ), d tSaP (0) G H S (R 3 ); (1.6) 

sup \g a p(0, x, x 2 , x 3 ) — m a p\ — > as r — > oo, (1.7) 

|x|=r 

/jere |x| = ^/xf + x| + x| and m Q/ g is i/ie Minkowski metric. 

Then for s > |, t/iere is a positive T > and a unique semi-metric g a/ g(i) which satisfies 
(1-4) an d such that 

(g af3 (t)-m a p)eC([0,T],H s+1 ) and d t g(t) e C([0,T], H s ). (1.8) 

In addition, if the pair (ga/?(0), dtg aj g(0)) satisfyies the constraint the equations (1.5), then 
the metric g a p(t) is a solution to the vacuum Einstein equation (1.1). 

S. Klainerman and I. Rodnianski [20] succeeded in improving regularity below the critical 
index |. 

Theorem B. Assume the conditions (1.6) and (1.7) of Theorem A hold and g a p{t) is a 
classical solution to (1.4) such that (g Q/ a(0), ^g a/ g(0)) satisfy the constraint the equations 
(1.5). Then for s > 1 i/iere is a positive T > depending on ||dg a/ g(0)||#s suc/i #ia£ 

II^WILfz- <C||9g^(0)||^. (1.9) 



1.2 Solutions of the constraint equations 

The spaces are inappropriate for solutions of the constraint equations, roughly speaking, 
because in these spaces the Laplacian is not invertible. It turns out that the Nirenberg- 
Walker-Cantor weighted Sobolev spaces H m ^ ([23], [5]) are suitable for asymptotically flat 
manifolds and indeed these spaces have been widely used in General Relativity. We denote 
the norm of the weighted Sobolev spaces by 

IMIL,. = E IK 1 + \^\) i5+H) d a u\\h m , -oo < 5 < oo, (1.10) 

|a|<m 

and the space H m $ is the completion of C^°(1R 3 ) under the norm (1.10). 
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We may express now the asymptotically flat condition for a Riemannian metric h ab by 
means of (h ab — e ab ) G H m $. In this paper the identity is denoted by e ab . 

There is an extensive literature for the solutions of the constraint equations in asymptoti- 
cally flat manifolds. The following theorem was proved under various assumptions in [1], 
[6],[10][8], [9], [12], [22]. 

Theorem C. Let m be an integer greater or equal to one, — | < 5 < — |. Given a 
set of free data (h ab ,K ab ) such that ((h ab - e ab ),K ab ) G H m+1;S x H m>s+1 . Then there 
are exists a conformally equivalent data (h ab , K a b) which satisfies the constraint equations 
(1.5). Moreover, there is a constant C such that 

II (h-ai — G ab , K ab ) || „ „ < C (h a 5 — e ab , K afe ) 



1.3 The main result 

The disadvantage of the present situation is the inconsistency of the Sobolev spaces of 
Theorems A and B with those of the constraint equations. The initial data for the semi- 
linear wave equations (1.4) are given in fP-spaces, while Theorem C provides the initial 
data (solutions to the constraint equations (1.5)) in H m s. Therefore it is impossible to 
obtain a solution to the Cauchy problem for the vacuum Einstein equations with initial 
data which are given in one type of Sobolev spaces. Our goal is to unify the Sobolev spaces 
of the constraint and the evolution equations. 

Before stating the main theorem we need to introduce the extension of the spaces H m> $ 
into fractional order. We denote a scaling with e by f e {x) = f(ex). 

Definition 1.1 (H s ^ Sobolev spaces) For s > and — oo < 5 < oo, we define the 
H S) s norm by 

3 

The sequence {ifjj} C C^°(IR 3 ) satisfies the following: ipj(x) = 1 on Kj = {x : 2 J ~ 3 < \x\ < 
2i+ 2 }, j = 1,2,..., K = {x : |x| < 4}; sup(^) C {x : 2^ 4 < |x| < for j > 1, 

sup(^o) C {x : \x\ < 2 3 }; \d a ijjj(x)\ < C a 2~^ , where the constant C a does not depend on 
j- 

The space H s ^ is the set of all temperate distributions having a finite norm given by (1.11). 
Triebel [26] proved that whenever s is equal to an integer m, then 

oo 

J2 2il+s)2j \\(^hA\ 2 H - - E IK 1 + M) (5+H) ^lli w - (i.i2) 

3=0 |o|<m 

Thus whenever the parameter s is an integer, the norms (1.10) and (1.11) are equivalent. 



4 



Theorem 1.2 (Main results) Let s > | and — | < S < — |. Given a set of free data 
(h ab ,K ab ) such that ((h ab - e ab ),K ab ) G H s+hS x if Sj(5+ i. 

(i) T/ien i/iere exzsfo a conformally equivalent data (h a b,K a b) which satisfies the con- 
straint equations (1.5). Moreover ((h ab — e a; ,),K ab ) G H s+ i t $ x H Si g+i and depend 
continuously on the norms of ((h ab — e ab ), K a &). 

(ii) T/ien i/iere exists a T > and a semi-metric g a p(t) solution to the vacuum Einstein 
equations (1.1) such that 



{Sap(t) - m a/3 ) G C([0,T},H s+1>s )nC\[0,T},H S7S+1 ) (1.13) 



and 



(ga/?(*) - m a 



■f»\\H, +li 



<C||(h o6 -e a6 ,K a6 )|L xH . (1.14) 



/or £ G [0,T]. TTie metric g aj g(£) i/ie unique solution to the reduce Einstein (1-4) 
with initial data (h a b,K a b). 

Remark 1.3 (Uniqueness) Since the Ricci tensor is invariant under diffeomorphisms, 
it is impossible to get a unique solution to the vacuum Einstein equation. Because if 
R-o/3(g) = and <p is a diffeomorphism, then the pull-back (p*ga/3 also satisfies (1.1). How- 
ever, it can be shown that if two metrics g a/ 3 andg a p satisfy (1.1) and (g a/3 — m a/3 ), (g a/3 — 
m «/j) £ H s+ i s s, then there is a coordinates transformation x a — > y a = f a (x^), which 
preserve the harmonic condition (1.3) and such that g^iv) = ga/3 ( x (y)) ff^ll^- This 
transformation can be established by means of solutions to linear wave equations (see [15], 
[10]) with coefficients in the H s ^-spaces. Thus we can apply the energy estimate Lemma 
4-3 and the tools of Section 2 to establish the existence of this transformation in H s+ \^. 
Previously this procedure has been applied with one more degree of differentiability, but 
recently, Planchon and Rodnianski found a trick which allows the obtaining of the diffeo- 
morphisms without losing regularity, see Section 4 in the monograph [13] for details. Thus 
we conclude that for asymptotically flats metrics which preserve the harmonic condition 
(1.3) the uniqueness holds up to a diffeomorphism. 

Remark 1.4 We would like to mention that the results of Christodoulou [11] and 
Christodoulou and O'Murchadha [12] differ from ours. They assume (h ab — e a b,K a b) 
G H s+1>s+ i x H s+lg+ 3, while the solutions (g a/ 3(t) — m a p) belong to H s+ i t s(Cle), where Qg 
is a certain unbounded region ofM. A . Thus in their setting, the rates of fall- off of the initial 
data and the semi-metric are different. In addition, they require the regularity condition 
s>3. 
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The idea to solve both the evolution and the constraint equations in the weighted Sobolev 
spaces of fractional order H s $ has previously appeared in [2] and [3] , but for the Einstein- 
Euler systems. The regularity condition for these systems is higher since they are coupled 
with a fluid. 

The plan of the paper is as follows. In Section 2 we present several properties of the 
fractional weighted Sobolev spaces. Section 3 deals with the reduction of the wave equations 
into a first order symmetric hyperbolic systems. The specific form of these hyperbolic 
system has an essential role in our approach. The energy estimates are established in 
Section 4. In Section 5 we treat the existence, uniqueness and continuity of semi-linear 
first order symmetric hyperbolic systems in the //^-spaces and the main result is proved 
in Section 6. In this paper Greek indices will take the values 0, 1, 2, 3 while Latin indices 
1,2,3. 

2 Weighted Sobelev spaces of fractional order 

Here we present the basic properties of these spaces and the equivalence between various 
norms. All these results were established in the appendices of [4] and [3]. At the end of 
the section we define a norm on product spaces. 

Definition 2.1 (Definitions of norms) 

• Let {ifij} be the sequence of functions in Definition 1.1. For any positive 7 we set 

3 

and we will use the convention \\u\\h sS1 = \\ u \\H aS - The subscripts 2 J mean a scaling 
with 2i , that is, (ip]u) 2 ]{x) = (^Ju) (2 J x) . 

• For a non-negative integer m and (3 6 K ; the space C™ is the set of all functions 
having continuous partial derivatives up to order m and such that the norm (2.2) is 
finite: 

\\u\\o* = SU P (C 1 + \x\f +]al \d a u(x)\) ■ (2-2) 

2.1 Some Properties of H s § 
Proposition 2.2 
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1. Equivalence of norms H s s and H s ^ a : For any positive 7, 

IMI 2 ^ = E 2(M5W IK^II^ - E 2(l+5)2i ll^»^H^ = IMIIw (2-3) 
3 3 

2. Equivalence of norms (2.1) and (1.10): For any nonnegative integer m, positive 
7 and 5 there holds 

= E 2(|+5)2J II^»^H^ - E + H) (5+W) ^«IIW (2.4) 

3 \a\<m 

3. H s 5-norm of a derivative: 

ll$ u lk- M+ i < Nk, a - (2.5) 

4. Algebra: If s±, s 2 > s, si + s 2 > s + | and Si + S 2 > 5 — §, i/ien 

HHk,* < Nk 3 ,, a - (2-6) 

5. Compact embedding: Lei < s' < s and 5' < 5, then the embedding 

H s ,s ^ Hs',5'- (2.7) 

compact. 

6. Embedding into the continuous: // s > | + m and 5 + | > t/ien 

INIc^ < C|Mk, 5 . (2.8) 

7. Third Moser's inequality: Let F : R m ^ R l be C h /+1 -function such that F(0) E 
H S) s and where N > [s] + 1 . Then there is a constant C such that for any u e H Sy s 

\\F(u)\\ Hs , < C\\F\\ CN+1 (1 + \\u\\ Hs , s + ||F(0)||^, (2.9) 

In particular, if s > §, then u 6 L°° and \\F(u)\\h sJ < C\\u\\h SiS + ||-F(0)|k >4 . 

8. Difference estimate: Suppose F is a C N+2 -function and u,v G H s s H L°° . Then 

\\F{u) - F{v)\\ HstS < C(\\u\\ L ~, \\v\\ L ~) (\\u\\ HaiS + \\v\\ H J \\u - v\\ HsS . (2.10) 

9. Density: 

(a) The class C^°(R 3 ) is dense in H Sy s- 
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(b) Given u G H S}S , s' > s > and 5' > 6. Then forp>0 there is u p G C~(R 3 ) 
and a positive constant C(p) such that 

\\u p - u\\ H , iS < p and \\u p \\ Hs , sl < C{p)\\u\\ HsS . (2.11) 

10. Mixed norm estimate: If u G H s ^ and d x u G H St s+i, then 

IMk +M ^ (IMIff.,* + Ru|k,, +1 ) • (2.12) 

The proof of (2.12) follows from the integral representation of the norm (1.11) (see [I,], [3] 
in the Appendix). 

The density property (b) where proved in [3], [4] for 5' = 5. Only a slight modification is 
needed to include it also for 5' > S and therefore we leave it to the reader. 

2.2 Product spaces 

Definition 2.3 (Product spaces) We set X S} $ = H s ^ x H St s+i x H St §+i, and the norm 
of a vector valued function V = (v\, V2, V3) G X St § is defined by 

IMI^HHI^ + Ik2||^ +1>2 + ||^ll i5+1 , 2 . (2.13) 

We will use the following convention: for a vector valued function u : R x M 3 — > M. N , we set 
U = (u, dtu, d x u), where d x u denotes the set of all partial derivatives in the space variable 
x G 1R 3 . Thus \\U\W =\\u\\% +\\d t u\\% + \\d x u\\ 2 H 
Essential of our approach is the following observation. 

Remark 2.4 It follows from the Mixed norm estimate (2.12) above that ifU(t,-) G X s ^, 
then 

Nv)lk +M <ll^(V)lk, 5 - 
3 First order hyperbolic symmetric systems 

The system of wave equations (1.4) can be transferred into a first order symmetric hyper- 
bolic system. The specific form of the hyperbolic system play en important role in our 
approach. 

Letting h aj g 7 = reduces the wave equations (1.4) into 

dtga/3 = h Q( ao 

d t h,so = z^u {2g 0a d a h, S0 + g ab d a h ySb + C^h^h^g^g^} (3.1) 

(_g00 ) -l gO 6^ h7<5a = (-g00)-l g a6 9ah7<50; 
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where the objects C^j^ are a combination of Kronecker deltas with integer coefficients. 
We would like now to write system (3.1) in a matrix form. We set 

where g, a/3 denotes the inverse matrix of g a/ g. By introducing the auxiliary vector valued 
functions 




U 



we can write the system (3.1) in the form 





a = 1, 2, 3, 



A°(u)d t U = ( Aa ( u ) + ca ) d aU + B(U)U, 



(3.2) 



a=l 



where 



A°(u) 



( 


eio 


Oio 


Oio 


Oio 







Oio 


eio 


Oio 


Oio 







Oio 


Oio 


u eio 


^ 12 eio 


9 13 




Oio 


Oio 


£ 21 eio 


~22 

9 eio 


g 23 


V 


Oio 


Oio 


£ 31 eio 


^ 32 eio 


g 33 



10 \ 

10 



(3.3) 



A a (u) 



and 



/ Oio 


Oio 




Oio 


Oio 




Oio 


\ 


Oio 


2g M e 10 


(g al - S al )e w 


(~g a2 - S a2 )e w 


(g a3 ~ ^ a3 )eio 




Oio 


(g al - S al )e 10 


Oio 


Oio 




Oio 




Oio 


(g a2 - S a2 )e w 


Oio 


Oio 




Oio 




V Oio 


(~g a3 - 5 a3 )e w 


Oio 


Oio 




Oio 


/ 






( Oio 


Oio Oio 


Oio Oio 












Oio 


Oio S al e w 


5 a2 e w 5 a3 e w 












Oio 


5 al ei Oio 


Oio Oio 












Oio 


5 a2 ei 10 


Oio Oio 












V Oio 


5 a3 ei 10 


Oio Oio 


J 







B(U) 



Here = (-g 00 )- 1 Cg^a e g c „g^g^. 



( Oio 


eio 


Oio 


Oio 


Oio 




Oio 






qkXI 






Oio 


Oio 


Oio 


Oio 


Oio 




Oio 


Oio 


Oio 


Oio 


Oio 




V Oio 


Oio 


Oio 


Oio 


Oio 


) 



(3.4) 



(3.5) 



(3.6) 
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Apart from the facts that A a {u) and C a are symmetric matrices and A°(u) is positive 
definite, they hold three additional properties: (i) the matrices A a (u) doest not depend on 
the derivatives of u; (ii) the coefficients of d t u in A°{u) do not depend on t; (iii) it follows 
from Moser type estimate 7 and Algebra 4 of Proposition 2.2 that if g a/ 3 — m a/ 3 G H s+ ij, 
then (g aa — 1) G i? s +i,<5 (a = 1, 2, 3) and g a/3 G H s+ ^s whenever a ^ (3. Thus the matrices 
(^4°(n) — e), ^4 a («) G H s+ i t s whenever g Q/ 3 — m a/ 3 G H s+ i t s, while C a are a constant matrices. 
These facts are crucial for the energy estimates. 



4 Energy Estimates 

We consider here energy estimates for a first order linear hyperbolic system of the form 
A°d t 



u 

d t u 
d x u 



3 

E 

a=l 



(A a + C a )d a d t u 



u \ 




1 « 


d t u 


M 


d t u 


d x u j 







(4-1] 



Here u : R x M 3 - ^ = (<,).. =123 , C a = «,.).. and B = (l>, ; ),, ,,, are 
5A^ x 5N block matrices having the sizes of their blocks according to the following structure 



N x N 


N x N 


N x 3N 


N x N 


N x N 


Nx3N 


3N x N 


3N x N 


3N x 3N 



(4.2) 



We assume the following conditions: 



a£ = for i ^ j; a° n = a° 2 



a 33 is symmetric and — v v < v a. 3S v < cqv v Vw G M ; 
A°(t, ■) - e G H S+1 , S ] 

A a are symmetric with = a" - = 0, a — 1, 2, 3; 
A a (t, ■) G H s+1 , s , a = 1,2, 3; 

C a are constant and symmetric with c" x = c" - = 0, a = 1, 2, 3; 
bj! = and bij are constant, i,j 
B(t, ■) := (bj J ) ij=2i 3 G 
.F(£, •) G H 8 ,s+i- 



1,2,3; 



(4.3a) 

(4.3b) 

(4.3c) 
(4.3d) 
(4.3e) 
(4.3f) 
(4.3g) 
(4.3h) 
(4.3i) 
(4-3j) 



Note that any system which is originated from a linearization of (3.1) meets the above 
requirements. 
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4.1 H s s - energy estimates 

We define an inner-product on X s ^ in accordance with the equations (4.1). Let 

A'(«) = (1 - A)f («) = T- 1 ((1 + |e| 2 )^) («), 
where J 7 denote the Fourier transform. 
Definition 4.1 (Inner-product) 

• Inner-product on L 2 : For vector valued functions f,g £ L 2 , we set 

(/,0)l* = / /V*, (4-4) 

where f T denotes the transpose matrix. 

• Inner-product on i/ Sj< 5: For Ui,0i £ -ff s ,«5 ; we set 



oo 

f3_ 



3=0 

Recall that the subscripts 2 J mean a scaling (see Definition 2.1). 

A weighted inner-product on H St s+i x H St s+i- For a matrix a° 3 which satisfies 
(4.3b) and (v 2 ,v 3 ), (0 2 ,0 3 ) e # s ,5+i x # s ,5+i> we set 

g 2 <i— ( A .(,?( : )) (;•) *-(^(i)r 



J=0 



Inner-product on X 8j s- For a matrix A which satisfies (4-3a)-(4-3b) and V = 
{vi, v 2 , v 3 ),$= (0i, 02, 3 ) e X S|5; we set 



M..A J.J r j./ s,o \ V ^3 / V 03 



s,5+l,a9„ 



We denote by ||V||x sa0 the norm which is associated with the inner-product (4.7). 
From condition (4.3b) we see that 
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\ v fx. SA „<c £ 2<f+i)2i H Ufa) 



1 1 2 
2i H^ s 



. j=0 



^ 2 (§+<5 + l)2j 
j=0 



(tfv 2 ) 2J \\ 2 HS + || {tfv 3 ) 2j 



C o{\H\h s , s>2 + IMIff., a+ i, a + INItf„ 4+1 , a } 



<*I|V||L,. 



Thus we have shown: 



Corollary 4.2 (Equivalence of norms) The norms which are defined by the inner- 
product (4-V and (2.13), satisfy 



V\\x.. < \W\ 



x 



S,A° 



< yfa\\V\ 



(4.9) 



For a vector valued function u(t, x), we set u(t) = u(t,x), U(t) = (u(t) , d t u(t) , d x u{t)) and 
the energy of U (t) is denoted by 

E(t) = ±(U(t),U(t)) XsSA0 . (4.10) 

The energy estimate in the product space X s> g is indispensable tool of our method. The 
next Lemma establishes it and its proof relies on tedious computations. The essential point 
is that fact that A a G H s+ i^- This enables to use the Kato-Ponce commutator estimate 
(Theorem 4.4) with the pseudodifferential operator A s d x rather then A s as in [4]. 

Lemma 4.3 (Energy estimates) Let s > |, S > — | and assume the coefficients of 
(4-1) satisfy conditions (4-3). IfU(t,-) G C^°(]R 3 ) is a solution to the linear system (4-1), 
then 

j t E(t)<Cc (E(t) + l), (4.11) 

where the constant C depends on \\(A° - e)\\ Hs+16 , ||-4°||^ S+M , \\B\\ HsS+1 , ||.F||ir M+lJ 
||5t-4.°||i°°> s and 5. 



An essential tool for deriving these estimates is the Kato & Ponce commutator estimate 



Theorem 4.4 (Kato and Ponce) Let P be a pseudodifferential operator in the class 
S? j0 , / G H s H C\ g G H s n L°° and s> 0. Then 

IIPC/^-ZP^IIy^cov/iUooyiH^ + ii/ii^yiioo}. (4.12) 
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Proof (of Lemma 4.3). Taking into account the structure of the inner-product (4.7) we 
see that 



1 oo 

The infinite sum of the right hand side of (4.13) is less than 



2J 



(4.13) 



L- 



oo 

f 3 



and 



^N\\dt (4s) \\L^Y. 2{ ~ 2+8+m W WM* II*- = V^II^II^H^II^ (4.14) 

3=0 

(u,d t u) SjS < \\u\\ Hs>s Jd t u\\ HsA2 < i (\\u\\ 2 HsA2 + ||ft«||lr., 4+lia ) • (4-15) 



We turn now to difficult task, this is the estimation of ( ( I , dt 



d t u \ g ( d t u 
d-u j ' 1 \ d x u 



s ,s+i,4 3 



Setting 

E 9t (j) = (a* (^d t u) 2J ,A s (^dt(d t u)) 2J ) L2 , (4.16) 

E 9x U) = (A s {^d x u) 2] , (a° 3 ) 2J A* (^dt(d x u)) 2J ) L2 , (4.17) 
and using the specific form of (4.6) we see that 

{ZM&^-t**™*™*™- (4 - 18) 

We first estimate Eg x (j). For that purpose we define a sequence of functions 

where {^j} is the sequence defined in Definition 1.1. This sequence has the following 
properties: ^ fe E C °°(M 3 ), \d a * k (x)\ < C a 2~\ EZo^k(x) = 1 and 

^ k (x)ipj(x)^0 only for k = j - 3, j + 4. (4.20) 
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We will use the convention that \1/ 



j-m 



whenever j — m < 0. Then 



E dm (j) = {A s [(^)J,(a° 3 ) 2j A 



,k=0 J 20 



2J 



L- 



j + 4: 



2J 



L 2 



(4.21) 



= £ {A s [^p x u) 2J \,{4 3 ) 2J A s [^ k i;p t (d x u)) 

k=j-3 
j+4 

k=j-3 



Our aim now is to take ag 3 across A s in (4.17), and then we can use the fact that U satisfies 
equation (4.1). In order to do it we will use the commutator estimate, Theorem 4.4. 
However, if use the commutator (4.12) directly with / = (^ k ) 2 j and g = {^ k ^dtd x u) 2:j , 
then we get || k ip?d t d x u) 2j \\ L «, < || (^ k ipjd t d x u) 2J \\ Hs by Sobolev inequality. That would 
leads to the condition s — 1 > |, and therefore we would not obtain the desired regularity. 
In order to avoid this, we write 



^ k tfd t d x u) 2J = -d x (^k^dtu) 23 - (d x (V k ip])) 2J (d t u) 2j , 



(4.22) 



then 



1 

2^ 
1 

2^ 
1 

2^ 
A s 



■ d x {^ k tfd t u) 



2> 



-A s 



(d x (y k rf)) 2J (d t u) 23 



(A s d x ) \(y k ijp t u)J - (A s ^) \{^dtu) 2] 



(4.23) 



(* fc ) 2J (A s d x ) (tfdtu) 
(d x (y k ^)) 2J (d t u) 2j 



2J 



Inserting the last three expressions in each term of Eg x (j, k) in the right hand side of (4.21) 
results in 



1 (A* [(^d x u) 2] ] , (4 3 ) 2] \(A°d x ) \{^d t u) 
- (* k ) 2 i(A s d x ) \{^d t u) 2] 



2J 



L- 



2-i 



- (A 6 



A s 



, (a° 3 ) 2J (A s d x ) [(^d t u) 

WMJ > ( a 3s) 2 , AS [fa (*^|)) 2J {dtu) 2j 



(4.24) 



L 2 



L 2 



--■ E d x (a, j, k) + E dx (b, j, k) + E dx (c, j, k) . 
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Estimation of Eg x (a, j, k): Applying the Cauchy Schwarz inequality we get 



\E dx (a,j,k)\ < 



2' 



' 23 1 1 L° 



X 



(A s d x ) {^ k tfdtu) 2J - (A s d x ) (tfdtu) 



2J 



L 2 



(4.25) 



The advantage of (4.22) is that A s d x e OPSffi, hence the Kato-Ponce (4.12) is applied 
with s + 1 rather than s. Therefore 



<l|V(**) a i|| L .||(^«) a ,|| jff . + ik^ii 



L 2 

H s + 1 23 



(4.26) 



From the properties of ipj (see Definition 1.1 ) and (4.19), we see that ||<9 a (\l/fc) 2 j < C, 
where the constant C is independent of j and k. Hence both || V(^ r fe) 2 i || ioo and || (^/c^j 
are bounded by a certain constant independent of j and k. For s > § the Sobolev inequality 



yields | (i/.'yC^ 



2J 



< 



|(^|<9 t -u) 2 j and combining these with inequality (4.25) we get 



\E dx (a,j,k)\< || (a°,) 



2J 1 1 L° 



'2.' 



(4.27) 



Estimation of E dx (c,j,k): Since (d x (^ k fpj)) ^ (d t u) 2 j = (F j:k ) 23 (ipjd t u) 2 j and the partial 
derivatives of {Fj )k ) 2 j up to order [s] are bounded by a constant C independent of j and k, 
there holds 



l^(c,i,fc)|<||(a° 3 ) 



2> 1 1 L° 



(iijd x u) 



2J 



H 3 



l jd t u] 



2J\\h° 



(4.28) 



Estimation of Eg x (b,j, k): We see from (4.24) that in order to use equation (4.1) we need 
to commute (^^33)21 with (A s d x ). Therefore we write 

= (A s d x ) [(rfd t u) 23 ] - (A s d x ) [(* fe a° 3 ) 2i (^ftu)*] 



+ (A s d x )[{* k 4 3 ) 2J (^d t uh] 

= (A s 4) - (A s d x ) [(Mis)* (tfdtu)*] 



(4.29) 



+2 j A s 



d x (y k 4^) 23 (d t u) 23 



($& 



2.' 



Thus (6, j, k) is sum of three terms: E dx (b, j, k) = E dx (d, j, k) + E dx (e, j, k) + E dx (/, j, k) . 
The first one will be estimated by means Theorem 4.4, in the second one we use algebra 
property of H s and in the last one brings us to equation (4.1). 
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L 2 



We recall that (A s d x ) 6 OPS^ 1 , therefore by Kato- Ponce commutator estimate (4.12) 
II ( A ^,) Mdtuh] - (A s d x ) {tfdtu)*] 

<{||V(^a° 3 ) 2J || ioo \\(^d t u)»\\ Hs + ||(^a° 3 ) 

From (4.19) and (4.20) we see that 



)23 Wffs+l \\{iiptu) 2 j\\ La 



(4.30) 



IVOMh) 



2J 



<Ci||(a° 3 ) 



2J 



C 2 2^||V(a° 3 ) 2J || LOO 



(4.31) 



and 



|(^a° 3 ) 2 ,||^ 



< C 3 ||(^a° 3 ) 2i || HS+1 = C 3 \\(ip k (a° 3 -e)) 2J + {^ k e) 2] \\ Hs+1 



= C 3 
<C 3 



((A (a° 3 -e)) 2fc ) 2J _ fc + ((^e) fc ) 
((^ (a°3-e)) 2fe ) 2J _ fc 



2.'- 



C 3 {||(^a° 3 -e) 



C 3 ||((^fce 



(4.32) 



kl2i~ k \\H S + 1 



2 k \\H S - 



< 



r oo rs_/ 



(4.33) 



Thus, the combination of (4.31) and (4.32) with the inequality || (■?/>! 
l^l^iM)^ !!^ (keeping in mind the factor 2~ j in (4.24)), leads to 

\E dx (dJ,k)\ 
= (A s , (^ fc a° 3 ) 23 (A s d x ) [(^d t u) 2J ] 

-{A s d x )[^A) 23 {^d t u) 2] ]) L2 

Z { 1 1 ( a °3) » ILoo + ll ( Va° 3 ) 2 J ioo + 1 1 (a° 3 - e)) 2fc 1 | } 

xiK^u^n^iK^ti^ii^.. 

We turn now to E dx (e,j,k) = (A s [{^d x u)^] , (A s d x ) [(* fe a° 3 ) 2 , (rfd t u) 2] ]) L2 . Noting 
that 

d x (^Atf) 2j (d tU )» 

=d x (* fc a 33 ^j) 2J {if}jd t u) 23 + 2 J (^f k al 3 d x ipj) 2:j (ipjd t u) 2j , 
and applying the algebra property of H s , we get 



(4.34) 



d x ^ k ^]) 2] {d t u) 2] 



< 



|^(^a° 3 ^) 2i || Hs ||(^u) 



2^|(^a° 3 ^) 2J || Hs W^AAAhs 



(4.35) 



Now lla^^a^),, L« < II {^k^i 3 ) 2 j \\ H s+l and I^Wfca^VjOaill^ £ 
hence by (4.35) and inequality (4.32) we get 



\E dx {e,j,k)\ 

ll(^(a° 3 - e)) 2k \\ Hs+1 + l} ||(^«)2*L. II(«^'IU • 



(4.36) 
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In order to use equation (4.1) we write 

\ 

A a 



/ 

V 



o 



A a 



C" 



( o 
V 











c a 



, a = 1,2, 3, 



(4.37) 



where A a = (a^)j J=2i 3, C a = (c^)j J=2i 3 are symmetric block matrix and c"- are constant. 
Further, let {^k} be the sequence which is defined by (4.19), then 



E dt (j) = (A s (ij]d t u)v,A s (i;p t (d t u)) 2] ) 



L 2 



A s (i(j d t u) 2J ,A* 



2J 



,fc=0 



2> 



L 2 



j+4 



2J 



L 2 



-- {^ s (^d t u) 23 ,A s [(Vktfdtidtu)) 

k=j-3 
j+4 

■ E E 8t (j,k). 

k=j-3 



From (4.21), (4.24) and (4.29) we see that 

E 9x (fJ,k) = (A s [(rfd x u) 2j ] ,A S [(^.a° 3 ) 2J ^d t d x u) v ]) 
and since U(t) satisfies (4.1), we have obtained 
{E 9t {j,k) + E ax (f,j,k)} = 



(4.38) 



L 2 



A" 



^ \d x u 



a=l x 

+ (a 



d t u 
I 



2J 



,A S 



e 

a° 3 



t 



d t u 
d r u 



2> 



L 2 



d t u 



d t u 
3 \d T u 



2J 



2> J 



+ A 



d t u 
d x u 



2J 



,A 



^k^[{A a + c a )d a r d 



d t u 



2 / u 22 "23 



Mb 



32 O33 



d t u 



.r 'I 



23 J / L 2 



2J J / L 2 



(4.39) 



2J 



L 2 



The main difficulty is the estimation of the first term of the right hand side of (4.39). We 
recall that C a are constant and A a G H s+ \s, therefore we may write 



* k tf (A a + C a )d a 



d t u 
d x u 



2J 



23 



jd a i*k^{A a +c a )r d 



d t u 



rU 



2> 



(4.40) 



2J 
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and hence 



A 
1 



2? 



2J 

i'd t u 



3 \d x u 



V k A a + C 



23 



(A s d a ] 



2 i d t u 



•p 



2J 



(4.41) 



d a U 



2J 



2^" 



The first term of the right hand side of (4.41) will be estimate by Theorem 4.4 with 
P = A s d a , in the second one the symmetry of A a will be exploited and in the third one 
we will use algebra property of H s . In both the first and third we take the advantage that 
A a e H S+1S . 



< 



+ 



(A s d a 

V(^ fe [A a + C a U23 
^ k (A a + C a )) 2J 



V k (A a + C a )) 2j (A s d a ) 



d t u 
d x u 



\\(ipp t u) 23 \\ Ha + \\(i/jp x u) 2j \\ Hs 

\($dtU)2j\\ TOO + \\(tfdxU)zi\\ TO , 



Now, for k = j — 3, + 4, 

V ( A a + C a 

and 



2.' 



< 



.4° 



+ 1 + 2 J ' 



2.' 



< 



(<M a 



'2-' 



2 k 



23- 



H s + 1 



1 < 



V4 a 



+ 1 



1 • 



In addition, since s > |, 



IWM2i|| LOO < IK^jW^H^. and IK^a-u)^! 



< 



I (tfd x u) 



2J 



L' 2 



(4.42) 



(4.43) 
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Thus, 



— I ( A s 
2i 







d t u 
3 \ d T u 



13 J 



, (A'0 a ) 



I a i na \ 

,2 I d t u 



V k [A a + C a V 



2-> 



V k [A a + C a )) 23 (A s d a ) 



d t u 
j \d r u 



2J 



I? 



< 



VA a 



A a 



(ip k A a ) 2 , 

2 



x \ \\(tfdtv>h\\ H . + \\(tppxu) 2 o u 



(4.44) 



As to the third term of (4.41), writing 



2 7a\ ( dtu 



d a (V k ^A a ) (q u ) = {daiVrts A a ) + 2(^ k A a d a ^) ^ U 



d t u 



and noting that 



2J 



< 



H 3 



2J 



H s+1 



< 



*PkA a 



2 k 



H s + 1 



then by the embedding H s+1 ^ H s we get that 



A s 



2 I d t u 
( d x u 
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A" 



< 



(AA a ) 



2 k 



+ 



H s+1 



2J 



L- 



x {\\(tfdtu) 2 j\\ Hs + 11(^1^)2,11^) {\\(tl}jdtu)2i\\ H s + \\(ipjd x u) 2 j\\ HS ) . 



(4.45) 



We turn now the second term of (4.41). Recall U is a C(f(IR 3 ), therefore A s (d t u), A s (d x u) 
are rapidly decreasing functions. This allows us to make the following operations: 



d a { A 



2 / d t u 
^ I d x u 



+ 



A s 



I 2 I d t u 



3 \d x u 



d t u 



21 



2J 



T 



^ k [A a +C a ) ) A 



(A a + C a ) ) A s 







d t u 
3 \d x u 

I d t u 



* \ flL„ 



h J < 
0. 



( d t u 
\d x u 



2J 



2J 



d a (^ k (A a +C a )) 23 A s 



3 \d x u 
I d t u 



* fc ( A a + C a 



A s 



2 > 



3 \d x u 



2 I d t u 



3 \d x u 



2J 



2J 



2J 



2J 



dx 



dx 



dx 



dx 



(4.46) 
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Since A a and C a are symmetric the first and the third terms of the right hand side of 
(4.46) are equal and hence 



2 

2^ 

1 

"2i 



v4 



d t u 
3 \d r u 



d t u 



3 \d x u 



2J 



*k [A a + C 



A s 



2J 



A* 



2J 



a / 2 ( dtu 



3 \d:rU 



2> 



21 



L- 



L' 1 



(4.47) 



< 



A a 



2J 



+ 



( ) 2i L°° 



That completes the estimation of the first term of the right hand side of (4.39). The second 
and the third terms of are easier to handle since they do not contain derivatives of high order. 
Recalling conditions (4.3i) and (4.3j) and using algebra in H s for s > |, we have 



< 



A' s 







d t u 



3 \ 8 



2J 



3 \ b 32 b 33 



d t u 



2J 



L 2 



2J 



2i 



+ 



(4>]d x u) 



2J 



(4.48) 



X (W^A^vWh* + IK^A^IIff*) 



and 



< 



A s 



d t u 
d x u 



2J 



,A S 



h 
h 



2J 



'2J 



L- 



2.' 



(4.49) 



To complete the proof we need to summarize 

oo / i+4 



£ £ 2(i^+^(E ft 0-,A ! )+^(7,fc)))- 

j=0 Vfc=J'-3 



(4.50) 



We see from the inequalities (4.27), (4.28), (4.33), (4.36), (4.44), (4.45), (4.47), (4.48) and (4.49) 
that the estimation of (4.50) consists of the following types of series: 



Type 1: 



oo / J+4 



j=0 \k=j-3 

where / and g belong to H s ^+i and h is in L°°; 



(4.51) 



Type 2: 



oo / j+4 



i=0 \k=j-3 



)2»IU» 11(^7)2^11^ II(^)2»IIh' 
where / and belong to H s 5+1 and /i is in L°°; 



(4.52) 
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Type 3: 

oo / j+4 \ 

52 52 ^ +&+m (\\{i>kh) 2k \\ H s + . Mf)A H s \M 9 )A H s) , (4-53) 

3=0 \k=j-3 J 

where / and g belong to H s ^+\ and h is in H s+ \^; 

Type 4: 

oo / J+4 \ 
52 E ^ +5+m (11(^)2*11^ Wif)A H ' W^9)»\\ H s) , (4-54) 

3=0 \k=j-3 J 
where / and g belong to H s s+i and h is in Hg+i^; 

Type 5: 

oo / i+4 \ 

52 E ^ +&+m (\\^jh)v\\ H s + i \Mf)A HS \Mg)A H ) > ( 4 - 55 ) 

i=0 \fc=j-3 / 

where / and g belong to H St s+i and h is in ffg+i^; 

Type 6: 

oo / i+4 \ 
£ 52 2 (f +5+1 ^ (iK^^H^i ||(^ 2 /) 2 ,||^ IIGfoM*.) , (4-56) 

3=0 \k=j-3 J 

where / G iJ Sj< 5, # G and /i G H s ,s+i- 

The estimation (4.51)- (4.56) will be done by means of the Cauchy-Schwarz and Holder's inequal- 
ities and the equivalence property (2.3), of the H Sj $-norm.. Starting with type 1, we see from 
(4.51) that 



oo j+4 

l3 



52 E 2 ( ^ + ^(ll^ 2 /) 2 4^IK^) 2 HI^ 

j=0 k=j-3 



< 7\\ h \\ L ~522(i^ (\Mfh4l + IK^MlU (457) 



j=0 

l^(ll/ll^ +1 , 2 + ll^llk W 



Similarly we estimate type 2, the only difference is the use of the equivalence (2.3) in the final step. 
Number 3 is more sophisticated, we first note that (§+<5+l)2j < ^-\-S)j+(^+5+l)j + (^+5+l)j 
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for 8 > — |, then we apply the Holder inequality with i, | and \ and get 



oo / J+4 



j=0 \fc=j-3 
oo j+4 



< E E 2(|+5)J IKV-^llfl^x) (2ti^ ik^/mu 2( 



+S+l)j 



2J 



j=Q k=j-3 



< 



oo j+A 

E E 

, i=0 fc=j-3 



oo j+4 

E E ( 2 

j=0 fc=j-3 



oo i+4 

E E ( 2 <« 

j=0 k=j-3 



3 +S+l)j " 



IK^H 



( ** +1 »"|Wfr)* 1L 



(oo i+4 
E E (2 { ^ +<5)2fc 
j=0 fe=i-3 



•h) 2 k \\ H s 



E(^ +s+1)2j \\m*\\l 



^(2(I + W||(^) 2 ,||^ 



2 k \\H a + l 



E( 2(|+s+w II 



k J=0 



t+^IK^) 



< 



H 



s+1,5 



H 



3,5+1,2 



\H. 



s,5+l,2 



< 



H 



s+1,5 
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s,a+i,2 



+ 



\H 



a,8+l,2 



The estimations of Types 4, 5 and 6 are similar to the last one. 
We may conclude now that 



d t u 
d x u 



dtu 
d x u 



oo / j+4 

£ E ^ +5+1)2j (E dt ti,k)+E dx (j,k)) 

j=0 \k=j-3 



< C 



+ 1 < C 



x +1 

-X-s,5 



(4.58) 



where the constant C depends on ||(-4° — e)||# s+1 s , \\A. a \\H s+1 ji \\B\\h SiS+1 , 

II^x^^IIl 00 ; s and 5. By the embeddings Proposition 2.2:6 and 2.2:3, we may replace ||.A a ||£oo and 

||3:r.4 a ||,L°° by their corresponding H s g norm. Thus combining inequalities (4.14), (4.15), (4.58) 
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with Corollary 4.2 we get that 



~{U(t), U(t)) Xs SA() < C (\\Uf Xstg + 1) < Cc (\\Uf XsAA0 + 1) , (4.59) 

here cq is the constant of the equivalence (4.3b) and in addition C also depends on ||9ta.3 3 ||ioo = 
||9^°||l°o. This completes the proof of the energy estimates. □ 

4.2 L 2 5 - energy estimates 

The L| space is the closure of all continuous functions with respect to the norm 



^ = I (l + \ x \y d \u(x)\ 2 dx. (4.60) 



Similarly to Definition 2.3, we set Yg = x L"j +l x L$ +l and the norm of V = (i>i, i^, 1*3) £ Y$ 
is denoted by 

nir 5 = lbi|l!| + ll^||| |+i + INHi, +i . (4.61) 
The equivalence of norms 1 1 V [ | 5 — ||V||y 4 follows from Proposition 2.2:2. 

In analogous to Definition 4.1, we define an inner-product which is appropriate to the system 
(4.1). So let ag 3 be a positive definite matrix and V, $ £ Y$ be two vector valued functions. We 
define an inner-product: 



(4.62) 



dx, 



and the norm which is associated with this product: = (V, V) T 2„o . If aSo satisfies 

(4.3b), then 

^ll^ll^, a o 3 <ll^||^<co||y||^ a o 3 . (4.63) 

Lemma 4.5 (Lg energy estimates) Assume the coefficients of (4-1) satisfy conditions (4-3a), 
(4-3b), (4-3e), and (4-3h). IfU(t,-) = (u(t,-),dtu(t,-),d x u(t,-)) £ X\ $ is a solution to the linear 
system (4-1), then 

J t \W)W 2 y s A 3 ^ Cc ° (ll^WHn,* + ' (464) 

where the constant C depends on the L 00 -norm of A 01 , d a A a and B. 
Proof (Lemma 4.5). Taking the derivative of (U(t), U(t)) L 2 a o yields, 

5 ' 33 

Ijt^K,^ = (U{t)AU{t)) YsAa + + \x\) 25+2 (d x ufd t 4 3 (d x u)dx. (4.65) 

By the Cauchy Schwarz inequality, the second term of the right hand side of (4.65) is less than 

Vm\\d t 4 3 \\ L ^\\d x u\\ 2 Ls+i . (4.66) 
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Let A a and C a be the matrices which is defined in (4.37), since U(t) satisfies system (4.1) we have 



(U(t),d t U(t)) Ys>a o a =1(1 + \x\) 25 u T (d t u)d 



a=l J 



25+2 



{{d t u) T , (d x uf) (A a + c a )d a 



d t u 

da-U 



dx 



+ /(i + N) 



25+2 



{(d t u) T ,(d x uf 
+ J(l + \x\) 25+2 [(dtufh + [d x u) T h] dx 

3 

= : Li + ^L 2 , a + L z + L A . 



( b 2 2 


b 23 \ 


f d t u V 


V b 32 


b33 / 


\ d x u J 



dx 



a=l 



We are estimating each term separately: 

\Li\ < \\u\\ Lj \\d t u\\ q < l - (\\uf q + \\dtuf Ll+ ) , 

|2L 2 , a | < 2v^iv(|25 + 2|p a || L o + \\d a A a \\ L - + l) (||d t u|| 2 2 + \\d x u\\l* 
|L 3 | < 4iV||B|| L =o (||^||^ +i + \\d x u\? L * +1 

and 



1^1 < 77 M2 ?ii + \\h\\U, + \\d x u\\l*, + 11/3112? 



5+1 



".5+1 



"5+1 



In (4.69) we have used the identity 
0= / 9a |(l + |x|) 25+2 
(2<5 + 2)(l + |s|) 



((d t u) T ,(d x u) T ) (A a + C a 



d t u 
d x u 



dx 



+ J(l + \ x \) 2 ^ 2 
+ j\l + \ x \) 25 + 2 
+ l\l + \ x \) 2S + 2 



25+1 x a 

X 



((d t u) T ,(d x u) T ) (A a + C 



d a ((dtu) 1 ,(d x uY )[A a + C a 



d t u 
d x u 



(d t u) T ,(d x uf) [A a + C a )d a ( 2 



d t u 



r u 



d t u 
d x u 

dx 
dx 



dx 



{d t uf ,(d x uf)d a A 



d t u 
d x u 



dx. 



and exploited the symmetry of A a and C a . 



(4.67) 



(4.68) 

(4.69) 
(4.70) 

(4.71) 



Summing the inequalities (4.66), (4.68), (4.69), (4.70) and (4.71) and taking into account the 
equivalence (4.63 ), we get inequality (4.64). □ 
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5 Local Existence of Quasi-linear Hyperbolic Systems 



Let « : K x I 3 -> R N and set U = (u, dtu, d x u), we consider a quasi-linear first order hyperbolic 
system 

3 

A°(u)d t U = ^ {A a {u) + C a ) d a U + B(U)U (5.1) 

a=l 

under the following conditions: 

Assumptions 5.1 All the matrices are smooth function of their arguments and 

1. A°(u) , A a (u) and C a are symmetric matrices; 

2. A°(u) = (8Lij(u))ij=i ; 2,3 is a block matrix such that a^-(u) = for i / j and a^(u) = e for 
i = 1,2; 

3. A a (u) = [a1j{u))ij = i t 2,z are block matrices such that af(u) = for j,a = 1,2,3; 

4- C a = (cfj)ij = i^,3 are constant block matrices such that cj- = for j,a = 1,2,3; 

5. B(U) = (bji(C/))j, = i 2 3 is a block matrix such that hn(U) = and bi j(U) are constant, 
i,j = 1,2,3. 

The sizes of the blocks are ruled according to (4-2). 

Clearly the system (3.1) satisfies these assumptions. The main result of this section is the well- 
posedness of the system (5.1) in X Si< 5-spaces. 

Theorem 5.2 (Well-posedness of quasi-linear hyperbolic symmetric systems) Let 

s > |, 5 > -§, (f,g) G H s+ i tS x H StS+ i and suppose 

-v T v < v T &l 3 (f)v < fiv T v, \/v£R 3N and some u G R + . (5.2) 

Then under Assumptions 5.1 there exits a positive T a unique U(t) = (u(t) , dtu(t) , d x u{t)) a 
solution to (5.1) such that U(0,x) = (f(x),g(x),d x f(x)) and 

UeC([0,T\,X StS ). (5.3) 

Remark 5.3 We may conclude by Mixed norm estimate 10 of Proposition 2.2 and (5.3) that 

u G C([0, T],H s+ltS ) n C\[0, T],H S) s+i). (5.4) 

We adopt Majda's method and construct the solution through an iteration procedure [21]. Similar 
approach was carry out in [4], [3] for s > |. Here we will examine how the special assumptions 
of (5.1) enable us to improve the regularity. 
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5.1 Construction of the iteration scheme 



We first note that the Embedding 6 of Proposition 2.2 implies that the initial data (/, g, d x f) are 
continuous, hence there is a constant Co > 1 and a bounded domain G2 C containing / such 
that 

— v T v < v T a^ 3 (u)v < c$v T v for u £ G<i- (5-5) 
Co 

According to the density properties of H s $ (Proposition 2.2:9), there are sequences 
{/ fc }fcLo' {9 k }T=o c an d a positive constant R such that 

||(/V,^/°)|U s+M < C\\(f,g,d x f)\\ Xs ,, (5.6) 
l|«-/°lk,, +1 , a <R^uEG 2 , (5.7) 



and 



" (f,9M\\x.j < 2~ h £-. (5.8) 



The iteration scheme is defined as follows: Let U°(t,x) = (f°,g°,d x f°) and U k+l (t, x) 
(u k+1 (t,x),dtu k+1 (t,x),d x u k+1 (t,x)) be a solution to the linear initial value problem 



A°{u k )d t U k+1 = £* =1 {A a + C a ) {u k )d x U k+1 + B{U k )U k+l 
C/ fc+1 (0,x) = (f k (x),g k (x),d x f k (x)) 



(5.9) 



The linear theory of first order symmetric hyperbolic systems (see e.g. [18]) guarantees the 
existence of a sequence {U k (t)} C Co°(R 3 ). Therefore for each k 

T k = su V {T: sup \\U k (t)-(f ,g ,d x f°)\\ XsS < R} > 0. (5.10) 

0<t<T 

We claim that there is T* > such that T k > T* for all k. 



5.2 Boundedness in the X s j-norm 

Lemma 5.4 (Boundedness in the norm) There is a positive constant T* such that 

sup{T : sup \\U k (t) - (f°,g ,d x f°)\\x sS <R}>T* for all k. (5.11) 

0<t<T 

Proof (of Lemma 5.4). Let V k+l = U k+l — U°, then it satisfies the linear system 

3 

A\u k )d t V k+l = (A a {u k ) + C a ) d a V k+1 + B(U k )V k+1 + T k , (5.12) 

a=l 

where 

3 

= Yl ( Aa ( uh ) + c °) d aU° + B{U k )U° 

a=l 
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and V k+1 (0,x) = (f k+1 (x),g k+1 (x),d x f k+1 (x)) - (f°(x),g°{x),d x f°(x)). At this stage we need 
to verify that the linear system (5.12) meets all the requirements of the energy estimates Lemma 
4.3. Clearly the matrices A°(u k ), A a {u k ) and B(U k ) satisfy conditions (4.3a), (4.3b), (4.3e) and 
(4.3h). 

We check now that rest of the conditions of (4.3). From the induction hypothesis (5.10), we 
have that \\u k — f®\\ 2 H S2 + \\d x u k — d x f°\\ 2 H M 2 < R 2 , therefore by Proposition 2.2:10, 



\u 



f°\\H s+1S 2 — CR. Applying the equivalence (2.3) and Moser type estimates (Proposition 2.2:7), 
we have 

PV) - e\\ Ha+1>s < CxII^Hh^ < C X (\\u k - f\\ Hs+1A2 + ||/°lk +1> , 
KC^CR+WfWH^). 

Similarly we get ||w4 a (u fc )||# s+1 s 1 < C2 (CR + ||/ ||_ff s+1 {2 )- Here the constants C\ and C2 depend 
on ||?x fc H^oo , and ||-4° — e\\c m (G 2 )' ll"4 a |lc m (G 2 ) respectively, which implies that conditions (4.3c) and 
(4.3f) hold. Having shown (4.3c) and (4.3f), we conclude from Proposition 2.2:6 that A a (u k ) € 
((3 > 0). Combing it with inequalities (5.7) and (5.10) we get 

\\d t A\u k )\\ L ~ <sup|^-(«)|||^ U fc |Uoo <C\\d t u k \\ HsS+12 

G 2 OU 

< C(\\d t u k - g°\\ Hs!W + \\g°\\H 8 , +1 J < C(R+\\g°\\ HsW ), 

this gives condition (4.3d). In order to verify condition (4.3i), we denote by B{U k ) the non- 
constant blocks of B(U k ). Then we apply again Moser type estimates 7 and Algebra 4 of Propo- 
sition 2.2, together with induction hypothesis (5.10) and the structure of the matrix B yield 

\\B(U k )\\ H ^ +1 < C,\\U k \\ Xs:6 < C 3 (R+ ||(/VA/°)lk.,) . 



The constant C3 depends on C m norm of B(U) taking in a bounded region of R 5Ar and || ?7" & Hz, 00 - 
Finally, the H s ^ estimates of A a (u k ) and B(U k ) with Proposition 2.2:4 provide an upper bound 
for ||J rfc || Sj( 5 + i. Thus we have verified all the conditions (4.3). 

We conclude that the constant C of the energy estimate (4.11) depends only on R and the initial 
data \\(f ,g ,d x f°)\\ XaS , hence 

d_ 

dt' 

and the constant C of (5.13) is independent of k. By Gronwall's inequality 



±<V*{t),V\t)) x < Cc»{ly k ^V k { t)) XsAA0 +l}, (5.13) 



\W K {t)\\x sAA0 < e° Cot {\\V k (0)\\ x ^ AO+ Cc t) . (5.14) 
Taking into account condition (5.8) and Corollary 4.2, we get from (5.14) that 
™p \\V k (t)\\ XsS < 

0<t<T 

e CcoT {4 (\\(f k ,g k ,d x f k ) - (f,g,d x f)\\ XsS + \\(f°,g°,d x f ) - (f,9,d x f)\\ XsS ) + Cc t} 
< e Cc T f& +CcoT \ < R 2 
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provided that T < T* := sup{t : e Ccot ( ^ + Cc i ) < R 2 }. □ 



Having shown the boundedness of {U k } we may conclude by the Compact embedding, Proposition 
2.2:5, that U k — > U in the X s ' >( y-norm for any s' < s and 5' < 5. By the Mixed norm estimate 
10, u k — ► u in H s > + i t gi and if we chose | < s' < s, — § < 5' < S, then the Embedding into the 
continuous 6 implies that 

u k {t) -»u(t) in C^M 3 ), 
a 4 u fc (t) -► Sfe«(t), d x u fc (i) -» in C(M 3 ). 

Therefore = (u(t) , dtu(t) , d x u{t)) is a solution to system (5.1) for < i < T*. 

5.3 Weak convergence 

Here we show the weak converges of {U k } in X Sj g. We chose the simplest inner-product on X s> g, 
that is, for V = (vi,v 2 ,v 3 ),$ = {<t>i,<hAz) G -X" S)(5 , we set 

00 

oo 

+ E 2(f+m)2j ( A * .A' (l&fa)*)^ (5.15) 

j=0 

oo 

3=0 

This definition coincides with (4.7) in the case where A is the identity matrix. 
Proposition 5.5 Given < s < L ^-, S < ^f-, V G X s / t g> and $ G Xgi^s", i/ien 

KV,*)^,!^ 11^11^11*11^. (5.16) 

The proof of Proposition 5.14 appears in [4], [3] with 5' = 5" = 5. Only a slight modification of 
this proof is needed in order to include it to (5. 16). Therefore we leave it to the reader. 

Lemma 5.6 (Weak Convergence) For any $ G X Si s, 

]im(u k (t),$) =(U(t),^) XaS (5.17) 

uniformly for < t < T*. Consequently 

\Mt)\\x„ t <VmM\\U k (t)\\x., t (5.18) 

and hence the solution U(t) of the initial value problem (5.1) belongs to C w ([0, T*], X s ^), where 
C w denotes the space of functions which are continuous in the weak topology. 
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Proof (of Lemma 5.6). We recall that \\U k (t) — U{t)\\n , *, —> for s' < s and 5' < 5. We can 



(5.20) 



pick now s" and 5" such that s < s", s < ^4^-, 5 < 5" and 5 < Given $ G X s> s and e > 0, 

we may find, by Proposition 2.2:9, $<= G X s » g» such that 

II* - <Mx„, < 2l and ll*e||x s „ >5 „<C( e )||$||^, (5.19) 
where is the constant of (5.11). Writing 

(u k (t) - u(t), $) Xst = (c/ fc (i) - cr(t), $ e ) Xs>a 

+ <£/ fe (i) - U(t), ($ - $ e ) i =: 4 + JJ, 

we have by Proposition 5.5 and (5.19) that 

\I k \<\\U k (t)-U(t)\\ Xslsl C(e)\\<l>\\ Xs ^0. 

As to the second term of (5.20), since \\U k (t) — U(t)\\x aS < 2i? by (5.11), we get from the 
Cauchy-Schwarz inequality and (5.19) that 

\n k \ < \\u k (t) - u(t)\\ XaiS ||$ - d> e |U M < f| = 6. 

Thus, 



hmsup|(C/ fe (t)-[/(t),$) Xs J<e 
fc 

and this completes the proof of Lemma 5.6. □ 

5.4 Uniqueness 

Lemma 5.7 (Uniqueness) Suppose U(t),V(t) G X s 5 are solutions to the first order sym- 
metric hyperbolic system (5.1) with initial data (f,g) which satisfy (5.2), then U(t) = V(t). 

Proof (of Lemma 5.7). Put W(t) = U(t) — V(t), then it satisfies the linear equation 

/ A°(u)d t W = Y?a=i(A a {u) + C a )d a W + B{U)W + F ( . 

\W(0,x) = ' [b - 2l) 

where 

3 

T = (A°(u) - A°(v)) d t V + " d a V + (B(U) - B{V)) V. (5.22) 

0=1 

Since U £ X s> s, A a (u), dpA a (u) and B(U) are bounded, we can apply Lemma 4.5 and obtain 

(See (4.62) for the definition of the norm ||W||?,2 t J- 
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We turn now to the estimation of ||-? r |lx, 6+1 in terms of the difference \\U — VWy^&^u)- I^ 0111 
the structure of the matrices A a (u) in Assumptions 5.1, we see that (.A (it) — A°(v)) dtV = 
(a§ 3 («) - a§a(i;)) d t d x v and 

d t v 



where 



(A a (u) - A a (v)) d a V = (A a (u) - A a (v)j d a \^^ 
af 2 (p) a| 3 (p) 



A a (p) 



a = 1,2,3. 



Our idea is to use inequality (2.6) with s = 0, si = 1 and s 2 = s — 1 and then to apply the 
Difference estimate 8 of Proposition 2.2. We also note that by Proposition 2.2:2, ||tt — v\\ 2 H ^ ~ 

1 1 -zx — v\\ 2 2 + 1 1 (it — u)||?2 — \\U — V\\y • These yield the following estimations: 

L s L s+i 6 



- a 33^)) d tdxvf L 2 ^ || (ajj 3 (u) - a§g(u)) df^HliL 

o — |- 1 u, o + 1 

II 2 



|| { a 33( u ; 

o+l 

<C||(a° 3 («)-a°3W)fc M l|5Ai;||^ +2 

< C(lkl|ff s+M , ||w||h-.+i, 4 ) ||« - v\\ 2 Hltt \\9tv\\ 2 HsjS+1 
<C(\\u\\ Hs+h5 ,\\v\\ Hs+1 , s )\\V\\ 2 x ^ \\U-Vf Ys . 



Similarly, 



(>(ii)-A»)d a 



dtv 
d x v 



T 2 

^6+1 



< C(\\u\\ Ha+hS ,\\v\\H 3+hS )\\u - vf His (||^|||r. i4+1 + \\d*v\?H. tS+1 

< C(\\u\\ Hs+hS , \\v\\ Hs+hS ) \\V\\ 2 XsS \\U - Vf Ys . 
Writing B = B(p, q, r), then by Assumptions 5.1:5 we have 

(B(U) - B{V)) V = d t {u - v) ■ V q BV + d x {u - v) ■ V r BV. 
Hence the simple weighted L 2 estimate gives 



\\d t (u - v)V q BV\\U < \\V q B\\ 2 LOB \\d t (u - v)\\U 

o + l 0+1 



2 



and 



\\d x (u-v)V r BV\\ 2 L 2 < IIV^II^H^^-i;)!!^ \\V\\ 2 loo . 

o+l o+l 

Thus, inequalities (5.24)-(5.27) with the equivalence (4.63) show that 

< c \\v\\ XsS \\u - v\\ Yt < c \\v\\ XsS \\u - v\\ Y ^ Uu) 

Inserting (5.28) in (5.23) and using Gronwall's inequality we get that 

\\W{t)\\ 2 Y <e°**||w(o)||* 

and since W(0) = 0, it implies that W(t) = 0. 



(5.24) 



(5.25) 



(5.26) 
(5.27) 
(5.28) 



□ 
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5.5 Continuation in the norm 



Lemma 5.8 (Continuation in the norm) Let U(t) be a solutions to the first order sym- 
metric hyperbolic system (5.1) with initial data (f,g) which satisfy (5.2), then (5.3) holds. 



Proof (of Lemma 5.8). Since X St s is a Hilbert space it suffices to show that 



limsup ||C/(i)|| x 



< \\U(0)\\ 



x 



-5,-4° (/) 



Having proved the uniqueness, we may assume that U is the limit of the iteration sequence U k . 
Furthermore, since u k (t) — ► u(t) uniformly in [0,T*] and the matrix A depends solely on u, we 
see from the inner product (4.7) that for a given e > there is a positive integer ko such that 



IVI 



x_. 



<(l + e)||y|| x 



k>k , V £ X a x. 



(5.29) 



,6,A°(u(t)) s,5,-A (ufc(t)) 

Using the fact that u(t, •) — ► /(•) uniformly as t — > 0, Lemmas (5.6) and (4.3), and (5.8) we get 
limsup ^(i)^ = limsup ||tf(t)|& 



< limsup liminf 



U k+l {t) 



x 



< (1 + ef limsup liminf U k+1 (t) 



S,A°(u(t)) 
2 



< (1 + e) 2 limsup I liminf e 



Cc t 



X 

k+1 



U k+L (0) 



S,A°(u k (t)) t 
2 



X s,S,A°(u k (0)) 



+ Cc t 



< 



(l + ef\\U(0)\\x 



This completes the proof of the Lemma and thereby of Theorem 5.2. 



□ 



6 Proof of the main result 



The solution of the constraint equations (1.5) in the weighted Sobolev spaces of fractional order 
H Sj s has been proved by Maxwell [22] for s > ^ and Brauer and Karp for > 1 [3] (see also [2]). 
Thus for a given set of free data (h a &, K a fc) such that [h a b — e a b, K a fc) G H s+ u x H Si s+i, there is 
conformally equivalent data (h a ^,K a fe) which satisfies the constraint equations (1.5). Moreover, 
there is a constant C such that 



(h ab - e ab ,K ab )\\ Hs+i gxHs S+i <C||(h o6 -e a6 ,K o6 )||^ +ij5X ^ 5+i 



(6.1) 



We apply now Theorem 5.2 to (g Q/ g — m^, dtga/3, d x g a p) with initial data (1.2) and where the 
pair (h a fe,K a fe) satisfies the constraint equations (1.5). Then g aj g(t) is the unique solution to the 
reduced Einstein equation (1.4) and (1.13) holds by Remark 5.3. Inequality (1.14) follows from 
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(6.1) since for t G [0,T] the bounds of \\g a p(t) ~ ta a p\\H a+1>s and \\dtg a p(t)\\ Hs>s+1 depend solely 
on the initial data (h a fe,K a h). 

In order to assure that g a /?(t) satisfies the vacuum Einstein equation (1.1) we need to establish 
the harmonic condition (1.3). Recalling that F 11 satisfies the linear wave equation 

g#d a dpF» - Tl^d v F» = (6.2) 

(see e.g. [10], [27]), it thus suffices to show that F»(0,x) = d t F^(0,x) = 0. Hence by the 
uniqueness of linear hyperbolic systems, it follows that F^ = 0. Note that g a/3 — e a/3 G H a +igi 
Ta/jg ^ ^ H 8 ,s+l and s > |, therefore these facts allow us to use known uniqueness results for 
linear hyperbolic symmetric systems with coefficients in H s [15], [19], or alternatively, we apply 
the L^-energy estimate Lemma 4.5, combined with GronwalFs inequality. We can now use the free 
data dtgoa to get the condition F M (0,x) = 0. Then exploiting the fact that (h a b,K a b) satisfies 
the constraint equations (1.5) leads to the second condition dtF^(0,x) = 0, see e.g. [1], [27]. 
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